Very Extended X-ray and Ha Emission in M 82: Implications for 

the Superwind Phenomenon 



> 



Matthew D. Lehnert^ 
Max-Plank-Institut fiir extraterrestrische Physik, Postfach 1603, D-85740 Garching, 

Germany 

Timothy M. Heckman 

Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 2121^ 



<: 

oo ' and 



Kimberly A. Weaver 



, NASA Goddard Space Fhght Center, Greenbelt, MD 20771 and Department of Physics 

^ ! and Astronomy, Johns Hopkins University, Baltimore, MD 21218 

O 

o 

> 



^Visiting observer at Kitt Peak National Observatory of the National Optical Astronomy Observatories, 
operated by AURA under contract with the National Science Foundation. 



-2- 



ABSTRACT 

We discuss the properties and implications of the 3.5 x 0.9 arc minute (3.7 
X 0.9 kpc) region of spatially-coincident X-ray and Ha emission about 11 arc 
minutes (11.6 kpc) to the north of the prototypical starburst/superwind galaxy 
M 82 previously discussed by Pevine fc Bally 1999| . The total Ha flux from 
this ridge of emission is 1.5 x 10^^^ ergs s^^ cm or about 0.3% of the total 
M 82 Ha flux. The implied Ha luminosity of this region is 2.4 x 10'^^ ergs s~^. 
Diffuse soft X-ray emission is seen over the same region by the ROSAT PSPC 
and HRI. The PSPC X-ray spectrum is fit by thermal plasma absorbed by only 
the Galactic foreground column density {Hh = 3.7 x 10^° cm~^) and having a 
temperature of kT=0.80±0.17 keV. The total unabsorbed flux from the ridge is 
1.4 X 10~^^ ergs cm~^ s^^ (~ 2.2 x 10^^ ergs s~^), comprising about 0.7% of the 
total X-ray emission from M 82. 

We evaluate the relationship of the X-ray/Ha ridge to the M 82 superwind. 
The Ha emission could be excited by ionizing radiation from the starburst that 
propagates into the galactic halo along the cavity carved by the superwind. How- 
ever, the main properties of the X-ray emission can all be explained as being due 
to shock-heating driven as the superwind encounters a massive ionized cloud in 
the halo of M 82 (possibly related to the tidal debris seen in HI in the inter- 
acting M 81/M 82/NGC 3077 system). This encounter drives a slow shock into 
the cloud, which contributes to the excitation of the observed Ha emission. At 
the same time, a fast bow-shock develops in the superwind just upstream of the 
cloud, and this produces the observed X-ray emission. This interpretation would 
imply that the superwind has an outflow speed of roughly 800 km s~^, consistent 
with indirect estimates based on its general X-ray properties and the kinematics 
of the inner kpc-scale region of Ha filaments. The alternative, in which a much 
faster and more tenuous wind drives a fast radiative shock into a cloud that 
then produces both the X-ray and Ha emission is ruled out by the long radiative 
cooling times and the relatively quiescent Ha kinematics in this region. 

We suggest that wind-cloud interactions may be an important mechanism for 
generating X-ray and optical line emission in the halos of starbursts. Such inter- 
actions can establish that the wind has propagated out to substantially greater 
radii than could otherwise be surmised. This has potentially interesting impli- 
cations for the fate of the outflowing metal enriched material, and bears on the 
role of superwinds in the metal enrichment and heating of galactic halos and 
the intergalactic medium. In particular, the gas in the M 82 ridge is roughly 
two orders-of-magnitude hotter than the minimum "escape temperature" at this 
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radius, so this gas will not be retained by M 82. 

Subject headings: Galaxies: active — galaxies: nuclei — galaxies: starburst — 
galaxies: ISM — optical: galaxies — individual: M 82 



-4- 



1. Introduction 



The history of the study of M 82 is a long and rich one. M 82 is the prototypical 
starburst galaxy and since it is one of the nearest (D = 3.63 Mpc - [b'reedman et al. 1994D and 
brightest members of that class, it has been studied at almost every wavelength ranging from 



the 7-rays through to low frequency radio waves (e.g., [Bhattacharya et al. 199^ ; [Bregman 



Schulman, & Tomisaka 1995|; |Courvoisier et al. 1990; 


Lynds & Sandage 1963 




Shopbell & 


Bland-Hawthorn 1998 


; Larkin et al. 1994; Petuchowski et al. 1994 


; Seaquist, Kerton, & 


Bell 1994; 


Hughes et al. 1994; 


Renter et al. 1994|; |Moran & Lehnert 1997; Ptak et al. 1997; 



Strickland, Ponman, fc Stevens 1997| ). It is perhaps impossible for any paper, let alone one 



this brief, to do justice to this rich history. We refer the reader to excellent exposes on M 82 
in, for example, [Bland fc I'uUy 1988| and [lelesco (1988)| . 



M 82 is also the prototype of the "superwind" phenomenon ( [Axon fc Taylor 1978 



McCarthy, Heckman, & van Breugel 1987|; |Bland & 'iully 1988; 


Heckman, Armus, & Miley 


1990 


; Bregman, Schulman, & Tomisaka 1995|; ptrickland, Ponman, & Stevens 1997|; 


Shopbell 



Bland-Hawthorn 1998| ). These flows are powered by the collective kinetic energy and 
momentum injected by stellar winds and supernovae in starburst galaxies. While recent 
studies have established the ubiquity of superwinds in local starburst galaxies (e.g., [Lehnert fc 
Heckman 1996[ ; [Dahlem, Weaver, fc Heckman 199"B[ ), there is still only a rough understanding 
of either their dynamics or the processes responsible for the observed optical and X-ray 
emission (e.g., Strickland 1998[ and references therein). Without such understanding, it is 
difficult to quantitatively assess the role that superwinds have played in the formation and 
evolution of galaxies ( [Kautt'mann fc Chariot 1998 ; ^omerville, Primack, fc Faber 1998 ) and 
in the chemical enrichment, heating, and possible magnetization of the inter-galactic medium 
( pibson, Loewenstein, fc Mushotzky 199^ ; Ponman, Cannon, fc Navarro 1998| ; [Kronberg 
[Lesch, fc Hopp 1999[) . 



In this paper we discuss the optical and X-ray properties of the spatially-correlated 
region of X-ray and Ha line emission recently reported by Pevine fc Bally 1999[ (DB hereafter) 
to be located roughly 11 kpc above the disk of M 82, far beyond the region where the 
superwind has been previously investigated. We are motivated by the important implications 
of understanding such a region for the superwind theory. If this region could be shown to 
be physically related to the superwind that M 82 is driving, then it implies that winds 
driven even by energetically-modest starbursts like M 82 are able to reach large distances 
into the halo, making them more likely to be able to escape the gravitational potential 
of the galaxy and thereby supply mass, metal, and energy to the inter-galactic medium 
(IGM). Such a finding would have implications for our understanding of the evolution of the 
IGM, the halos of galaxies, and quasar absorption line systems. Our study of this emission 
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region also provides some unique clues about the physics that underlies the observed X-ray 
and optical emission from superwinds, and strongly suggests that wind-cloud collisions in 
starburst galaxy halos are an important emission process. 



2. Observations and Data Reduction 
2.1. X-ray Data 

The ROSAT PSPC data were obtained during the period from March through October 
1991 and had a total integration time of about 25 ksec. M 82 was placed approximately in 
the center of the field for all the observations, and the source size is small enough so that 
there is no obscuration of the X-ray emission from M 82 due to the mirror ribs. 

We extracted a spectrum from a region of diffuse X-ray emission located about 11 arc 
minutes north of M 82, and coincident with the region of Ha emission reported by DB. We 
will subsequently refer to this region as the "ridge" of emission. This area of X-ray and 
optical emission is oriented roughly parallel to the major axis of M 82, and was denoted as 



sources 8 and 9 in [Dahlem, Weaver, fc Heckman 1998| . We extracted source counts in an 



elliptical region with a major (minor) axis of 3.3 (1.8) arc minutes oriented in a position 
angle of 43°. This region is just sufficient to encompass all of the obvious X-ray emission, 
but excludes the soft point source at the ENE end of the ridge of emission that is discussed 
in §3.2 below. We have also extracted source counts for the portion of the ridge that is 
contaminated by this point source using a circular aperture with a diameter of 1.75 arcmin. 
The background for these spectra was taken from a large region where there were no apparent 
background sources to the south and west of M 82. The relative distance of the background 
region from the center of the PSPC field was approximately the same as that for the ridge 
of emission. The spectrum was then grouped such that each channel contained at least 25 
counts. The extracted background subtracted data were subsequently fitted with various 
models using the package XSPEC. 

We have also inspected the archival ROSAT HRI data for M 82. These data have lower 
sensitivity than the PSPC data, and provide no spectral information. However, they do 
provide significantly higher spatial resolution (~6 arcsec FWHM on-axis and ~15 arcsec at 
the position of the X-ray ridge). The data we have used comprise a 53 ksec exposure taken 
in April and May 1995. We have used these data only to study the morphology of the X-ray 
ridge. 



2.2. Optical Imaging 



Our optical images of M 82 were taken on the night of February 27, 1995 (UT) using 
the KPNO 0.9m telescope. We observed with a Tek 2048^ CCD and through the optics of 
0.9m, this yielded a scale of 0.68 arc seconds pixel"^ and a total field of view of over 23 
arc minutes. The night was not photometric, with cirrus around all night and intermittent 
patches of thicker clouds. M 82 was observed through two filters, a narrow-band filter with 
a full width at half maximum of 29A and a central wavelength of 6562A and a wider filter 
with a FWHM of 85A and a central wavelength of 6658A. The first filter effectively isolated 
Ha and underlying continuum and did not include emission from either [Nil] lines at 6548 
and 6584A. The second filter provided a line-free measure of the continuum and was used 
for continuum subtraction of the first narrow-band filter to produce a Ha line-only image of 
M 82. The total integration time in each filter was 1200 seconds. 

The images were reduced in the usual way (bias subtracted, fiat-fielded, and then the 
narrow-band continuum image was aligned and scaled to subtract the continuum from the 
narrow-band Ha image) using the reduction package IRAF01RAF is distributed by the Na- 
tional Optical Astronomy Observatories, which are operated by the Association of Universi- 
ties for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 

Since the data were taken under non-photometric conditions, fiuxing the data directly 
using observations of spectrophotometric standards is impossible. However, obtaining fiuxes 
and hence luminosities is very important to our analysis. Therefore, to estimate the Ha 
fiux and continuum fiux-density of the off-band filter of the continuum and emission-line 
nebulae of M 82, we bootstrapped from the fiux given in a 5.8 arc second circular aperture 
centered on knot C from |0'Connell fc Mangano (1978 )| . Knot C was chosen because it was 
easily identifiable in our images and it appears in a relatively uncomplicated region of line 
emission. This procedure yields a Ha fiux from M 82 of about 4.6 x 10~^^ ergs cm~^ s~^, 
which compares very well with the Ha fiux from M 82 estimated by [McCarthy, Heckman, fc| 
van Breugel 1987| of 4.5x10"^^ ergs cm~^ s~^ over a similar region 90" x 90". All subsequent 



results for the Ha fiux and luminosity of various features within the Ha image of M 82 will 
rely on this bootstrapping procedure. 

To estimate the continuum fiux-density we used the fiux-density of O'Connell & Mangano 
taken at 6560A. Since the center of our narrow-band continuum emission is 6658A, we are 
thus assuming that the continuum has a constant fiux-density between 6560A and 6658A. 
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To see if this extrapolation procedure is accurate, we compared the flux in R of M 82 in a 35 
arc second aperture of 1.97x 10^^^ ergs cm^^ s^^A"^ (from NED based on data originally 
obtained by [Johnson 1966| ). If we also take the total counts in the 6658A filter and determine 
the fiux-density by extrapolating the conversion factor from measuring knot C of O'Connell 
& Mangano, we find a fiux-density of 2.11 x 10~^^ ergs cm~^ s~^ in a 35 arc second 
aperture. Our extrapolated flux is within 10% (0.1 magnitudes) and thus this extrapolation 
is reasonably accurate. 



3. Results 

3.1. Optical Imaging 

The Ha and narrow-band continuum images of M 82 are quite spectacular (Figures 1 
and 2). We find that the line emission is very extended, with discernible Ha line emission 
extending to about 11 kpc above the plane of M 82. This very faint emission can be seen even 
more clearly in the deeper Ha image obtained by DB. Like previous authors (e.g., [Bland fc 



I'ully 19S^ ; [McCarthy, Heckman, fc van Breugel 1987[ ) , we find that the Ha line emission has 



a very complex morphology that is very suggestive of outflowing gas, with loops, tendrils, 
and filaments of line emission pointing out of the plane of the galaxy. 

An especially interesting feature reported by DB is the faint ridge of emission about 11 
arc minutes north of M 82 at a position angle of about —20° from the nucleus. This ridge 
is then about 11.6 kpc (in projection) north of M 82 at the adopted distance of 3.63 Mpc 
( [breedman et al. 1994[ ). The dimensions of the ridge down to surface brightnesses of about 
3.5 X 10~^^ ergs s~^ cm ~^ arcsec"^ are about 3.5 arc minutes long x 0.9 arc minutes wide, 
with the long axis being approximately parallel to the major axis of M 82. At the assumed 
distance of M 82, this size translates into approximately 3.7 kpc x 0.9 kpc and the total 
Ha flux from this ridge of emission is 1.5 x 10^^^ ergs s~^ cm This flux is reasonably 
close to the estimate in DB of 1.2 x 10^^^ ergs s~^ cm~^, and yields an Ha luminosity of 
2.4 X 10^^ ergs s~^. The total Ha flux for the whole of M 82 is approximately 4.6 x 10"^^ 
ergs s~^ cm Thus, the ridge of emission constitutes approximately 0.3% of the total Ha 
emission from M 82 (not corrected for internal extinction). While the ridge appears rather 
structureless in our Ha image, the deeper image published by DB shows a collection of bright 
knots and loops or filaments connected to and immersed in the more diffuse and pervasive 
emission. 



The narrow-band line-free continuum image reveals a complex morphology that is as 
interesting as the Ha image. Besides the obvious complex dust absorption that is seen in 
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the image (with some dark features oriented perpendicular to the plane of the galaxy), we 
also see filaments of continuum emission extended out of the plane of the galaxy (although 
the filaments are not entirely obvious in Figure 2). These filaments have twists and turns 
as seen projected on the sky, but generally point perpendicular to the plane of the galaxy. 
They extend outwards almost 3 arc minutes to the northwest and almost 2.5 arc minutes 
to the southeast - corresponding to a projected distance out of the plane of about 3 kpc. 
Even though the continuum emission is in many ways as complicated as that seen in the Ha 
image, we do not observe any continuum emission from the ridge of spatially coincident Ha 
- X-ray emission. Using the counts to flux-density conversion for the narrow-band line-free 
continuum image and the same aperture as used to measure the Ha flux (but with the point 
sources removed from the area), we find an upper limit for the possible continuum emission 
from the ridge of fainter than 15.8 R magnitudes. This then implies an absolute R magnitude 
fainter than —12.0. 

While we find no evidence for any diffuse, spatially-extended optical continuum associ- 
ated with the Ha ridge, we do find a relatively bright point source in the line-free image that 
is coincident with the X-ray point source seen at the ENE end of the ridge in the ROSAT 
HRl data. This point source has an R magnitude of 12.2. We will discuss its likely nature 
in the next section. 



3.2. X-ray Imaging and Spectroscopy 

The optical images described above were aligned with the X-ray images starting with 
the coordinates given for the X-ray data from ROSAT and using the positions of stars within 
the frame from the Guide Star Catalogue and using coordinates given in the header from the 
digital sky survey. We then slightly adjusted the X-ray pointing by identifying and aligning 
X-ray sources on the digital sky survey image of the region. The close spatial correspondence 
between the ridge of optical and X-ray emission about 11 arcmin to the North of M 82 is 
obvious (Figure 1). 

We do not display the ROSAT HRI image of the ridge, since these data show only 
that the ridge is largely diffuse and structureless with the 15 arcsec (260 pc) off-axis spatial 
resolution of the HRI. The only substructure is a point source at the easternmost end of the 
ridge. This source appears to coincide with a stellar object visible in Figure 2 (the relatively 
bright star-hke object located about 40 arcsec ENE of knot T in Figure 2 of DB). This is 
most likely a foreground star, as we will briefiy discuss below. 

Examination of the 0.25 keV, 0.75 keV, and 1.5 keV PSPC maps of M 82 presented by 
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Dahlem, Weaver, fc Heckman 1998| - hereafter DWH - shows that the X-ray morphology 



of the ridge is energy- dependent. The ENE portion of the ridge (source 9 in DWH) is 
relatively much more conspicuous in the 0.25 keV map compared to the rest of the ridge 
(source 8 in DWH). This reflects the contamination of the eastern-most ridge emission by 
the (presumably) unrelated point-source visible in the HRI image. 

As described in §2.1 above, we have therefore extracted a PSPC spectrum of the main 
body of the ridge from a region that excluded the point source. The relatively small number 
of detected photons precludes fitting complex multi-component models, and we fit a model 
of a single thermal ("MEKAL") component to the spectrum. We are motivated to use 
a thermal model by the detailed X-ray spectroscopy of the brighter portions of the M 82 
outflow (DWH; |lVloran fc Lehnert 1997| ; [Strickland, Fonman, fc Stevens The fits 



strongly preferred a low absorbing column, so we simply froze the column at the Galactic 
value of Nh = 3.7 x 10^° cm~^. The data quality were not adequate to fit for the metal 
abundance, so we froze this at the solar value. Changing the assumed metal abundance did 
not significantly affect the best-fit temperature, but does directly affect the normalization 
of the fit. The best-fit temperature was kT=0.76±0.13 keV (90% confidence interval with 
a reduced ~ l-Sj Figure 3). The relatively high temperature is similar to that measured 
in the ROSAT PSPC spectroscopy of the innermost region of the M 82 outflow (~0.7 keV; 
DWH; Strickland, Ponman, fc Stevens 1997). The total unabsorbed flux in the 0.1 - 2.4 keV 



band from this part of the ridge of emission is 8.5 x 10 ergs cm s . 

For the ENE part of the ridge - whose emission is contaminated by the point source 
described above - we have fit a two-component MEKAL model (Figure 3). Both components 
were frozen at solar abundances, one component (the ridge emission) had an absorbing 
column frozen at Nfj = 3.7 x 10^° cm~^, and the other (the stellar point source) had an 
absorbing column frozen at 0. This fit yields kT=0.85±0.17 keV for the ridge and kT~0.1 
keV for the point source. The corresponding unabsorbed fluxes are 5.8 x 10^^^ ergs cm~^ 
s~^ and 3.5 x 10~^^ ergs cm~^ s~^ for the ridge and point source respectively. Combining 
the fits for the two regions of the ridge, the total X-ray flux is then 1.4 x 10~^^ ergs cm~^ 
s~^. This comprises about 0.7% of the total unabsorbed X-ray emission from M 82 in the 
PSPC band (e.g., [Moran fc Lehnert 1997| ; DWH). The total X-ray luminosity of the ridge of 



emission is 2.2x10 ergs s (in the 0.1 - 2.4 keV band). 

As stated above, the ENE point source is identified with a stellar object having an R 
magnitude of 12.2 (AF^ = 1.6 x 10"^'' erg cm~^ s~^)- The X-ray/optical flux ratio suggests 
that the most likely nature of this source is a late-type dwarf star, and this would also be 
consistent with its soft ROSAT PSPC spectrum (e.g., plosner, Golub, fc Vaiana 1985|) . 
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4. Discussion 

4.1. Is the Ridge of Emission Associated with M 82? 

Being able to isolate the ridge of emission to the north of M 82 allows us to investigate 
its nature separately from the rest of M 82. Could this ridge of emission be due to another 
galaxy near M 82, or is it instead gas excited by the M 82 starburst? 

We have found that the ridge of X-ray and Ha emission has fairly high X-ray and Ha 
luminosities - both luminosities are in the range observed for star-forming dwarf galaxies 
(e.g., [Hunter, Hawley, fc Gallagher 1995| ; [Heckman et al. 1995| ). However, we see no evidence 



for any optical continuum emission from this ridge. The limit for the level of continuum 
emission implies that any putative galaxy must be very low luminosity, comparable to that 
of the local dwarf spheroidal galaxies ( [Mateo 1998 ). Moreover, the limit on the continuum 



emission implies an Ha equivalent width of greater than 180A. Examining the sample of 
galaxies with Mb fainter than —15 in the the Ha survey of nearby dwarf irregular galaxies 
by [Hunter, Hawley, fc Gallagher 199"3[, we find that these have typical global Ha equivalent 



widths of a few to a few tens of A. The only such galaxy with an Ha equivalent width 
comparable to our lower limit for the M 82 ridge is DD053. This galaxy has a ratio of HI 
mass to Ha luminosity (in solar units) of 420. The corresponding HI mass for the M 82 ridge 
would have to be 2.5 xlO^ Mq. The HI map in [Yun et al. 199^ implies a maximum HI mass 



in the M 82 ridge that is only about 10^ Mq. Moreover, while the M 82 ridge has an X-ray 
luminosity similar to that of the dwarf starburst galaxy NGC1569 ( [Heckman et al. 1995[ ), 



the M 82 ridge would have a ratio of X-ray to optical continuum luminosity that would be 
over 40 times higher than NGC1569. Thus, it is clear that the ridge in M 82 has properties 
that are very different from those of star-forming dwarf galaxies. This argument, together 
with the strong morphological connection between the M 82 ridge and the M 82 outfiow (see 
DB) , leads us to conclude that this ridge is not a dwarf star- forming galaxy. 

Other possibilities for this region, like a cloud in the halo of the Milky Way are also 
highly unlikely. The HI detected in the data of [Vun et al. 1993| and the Ha emission from the 
spectrum of DB have velocities such that it they are both likely to be associated with M 82 
and not the galaxy (Purton fc Hartmann 199^ ; Paffher, Reynolds, fc Tufte 1998 ; Reynolds 



[et al. 1995| ). Moreover, while galactic high velocity HI clouds often show soft X-ray excesses 
( [Kerp et al. 1996[ ; p<.erp et al. 1999| ), the associated X-ray emission is very soft (?»0. 1-0.2 



keV; see discussion in [Kerp et al. 199(j|) . 
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4.2. Photoionization of the Ridge by M 82 

The ridge is evidently related to, and excited by M 82. In the case of the Ha emission two 
possibilities are obvious. Either the superwind in M 82 has worked its way out to ~11 kpc 
above the plane of the galaxy and is shocking halo material that appears in plentiful supply 
around M 82 ([Yun et al. 1993| , 1994), or Lyman continuum radiation from the starburst is 
ionizing gas at very large distances from the disk. Only the first possibility could apply to 
the X-ray emission from the ridge. 

Let us suppose that ambient gas is being photoionized by the stellar population of M 82. 
The rate of ionizing photons from M 82 is approximately 10^^ s~^ ( [McLeod et al. 1993| and 



references therein). Using the projected size and distance from M 82 estimated for the ridge 
of emission given in the previous section, and assuming that the ridge is really a "box" 
with a depth equal to its projected length (i.e., 3.7 kpc x 3.7 kpc) we find that the ridge 
will intercept about 8 x 10^^ ionizing photons s~^. This of course is an upper limit since 
some (most?) of the ionizing photons will be absorbed or dust-scattered before they reach 
distances as great as 11.6 kpc out of the plane of M 82. We can estimate the rate of photons 
necessary to give the Ha luminosity of the ridge (2.4 x 10^* ergs s~^) by assuming that the 
gas is in photo-ionization equilibrium, i.e., the number of recombinations equals the number 
of absorptions. Using the numbers for the Case B effective recombination rate of Ha from 
(Jsterbrock 1989| we find that it is necessary to have a rate of ~ 2x 10^° ionizing photons 



s~^. This is ~40 times smaller than the upper limit for the number of photons that M 82 
is able to provide, so only about 2 to 3% of the ionizing radiation initially emitted in the 
direction of the ridge needs to make it to this distance. Such an 'escape fraction' is consistent 
with HUT spectra of starbursts below the Lyman edge ( [Leitherer et al. 1995| ; |Hurwitz~et 



al. 1997|) , and it seems plausible that the optimal path along which ionizing radiation could 



escape from a starburst is along the cavity carved-out by the superwind. 

DB have made similar arguments in favor of a starburst-photoionization origin for the 
Ha emission. However, photoionization by stars does not produce the hot X-ray emission 
that is coincident with the Ha emission. Instead, shock-heating by the superwind may 
be an efficient mechanism for producing coincident X-ray and line-emitting plasma. We 
now explore the energetic plausibility of this model, and revisit the issue of the role of 
photoionization in §4.5. 
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4.3. Shock Energetics and the Basic Physical Properties of the Ridge 

The ridge of emission has an X-ray luminosity of 2.2 x 10^^ ergs s~^ and a temperature of 
9x 10^ K. Assuming that the shock is strong and adiabatic, we have that Tpost-shock = Jq^i 
where /i is the mass per particle, Vg is the shock velocity, and k is Boltzmann's constant. 
Thus to produce gas with T = 9 x 10^ K requires a shock speed of 820 km s~^. Such a 
high shock speed is quite reasonable considering that M 82 has been estimated to have a 
deprojected outflow velocity in optical line emitting filaments observed along the minor axis 
of roughly 660 km s"^ QShopbell fc Bland-Hawthorn 199^ ), and the outflowing X-ray gas 
could in principle have a velocity as high as 3000 km s"^ ( [Chevalier fc Clegg 1985|) . It would 



take material traveling at 820 km s ^, ;^1.4 x 10^ yrs to travel the 11.6 kpc from M 82 to 
the ridge of emission. Since this is approximately the estimated lifetime of the starburst in 
M 82 ( [Rieke et al. 1^^ ; Pernloehr 199^ ; [Schreiber 19981 [Satyapal et al~Tgg7| ), it is very 



plausible that the starburst could drive a shock over its lifetime this far out of the galactic 
plane. 

If we assume that we are observing a "box" of emission with dimensions of 3.7 kpc x 3.7 
kpc X 0.9 kpc, the total volume is 3.7 x 10®^ cm~^. Assuming this region is occupied by hot 
gas flUing a fraction ix of this volume, the measured luminosity in the ROSAT band pass of 
0.1 to 2.4 keV implies that the X-ray emitting gas has an electron density, Ug^x = 5.4 x 10^'^ 
fx"^"^^ cm~^, pressure, Px ~ 2ne^xkT=1.3 x 10^^^ fx^^^^ dynes cm~^, cooling time, tcooi ~ 

yrs, a mass, Mx=2.0 X 10^ fx^/^ M0, and a total thermal 
energy of 8x 10^"^ fx^''^ ergs. Similarly, assuming that conditions of Case B recombination 
and a temperature of 10^ K holds in the region producing the Ha emission, and that this 
region also occupies a volume 3.7 x 10^^ cm~^, then the implied density and mass of this 
gas is He^Ha = 4.3 X 10~^ f^l"^ cmr'^ and M^q, =8.0 x 10^ f^^ Mq respectively. If we make 
the physically reasonable assumption that this cooler gas is in rough pressure balance with 
the hot X-ray gas, then the minimum pressure of 1.3 x 10~^^ dynes cm~^ would require that 
the Ha emitting gas is highly clumped {iua ~ 8 x 10^^). This in turn implies that nf. Ha = 
5 cm~^ and M-Ha =1-1 x 10^ Mq. We note that a small volume fllling-factor for this gas is 
consistent with the complex knotty, fllamentary structure seen in the deep Ha image in DB. 



The total bolometric luminosity and ionizing luminosity of M 82 ([McLeodetal. 1993| and 



references therein) can be used together with the starburst models of [Leitherer fc Heckman 



1995)1 to predict the rate at which the starburst injects mechanical energy and momentum. 
While these predictions will depend to some degree on the evolutionary state of the star- 
burst (cf. Figures 66 through 69 in Leitherer & Heckman) we estimate that the mechanical 
luminosity and momentum flux provided by M 82 are ~ 2.5 x 10^^ erg s~^ and ~ 2 x 10^^ 
dynes respectively. For this wind mechanical luminosity, the rate at which the X-ray ridge 
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intercepts this energy is 2.5 x 10^^ ^w~^ erg s"^, where fl^ is the total sohd angle into 
which the wind flows < An). Comparing this heating rate to the X-ray luminosity of 
the ridge (2.2 x 10^^ erg s~^) shows that the wind could easily power the emission. It could 
also supply the ridge's estimated thermal energy in a timescale of 10*^ fx^^^ years (much 
less than the wind/starburst lifetime). 

Matching the observed high pressure in the X-ray ridge is more challenging. A wind 
momentum flux of 2.5 x lO^'* dynes leads to a wind ram pressure at a radius of 11.6 kpc 
of 2 X 10~^^ ^w'^ dyne cm~^. Consistency with the X-ray pressure estimated above would 
require that the X-ray gas has a large volume fllling factor {fx ~ unity) and that the M 82 
wind is rather well coUimated {Qy^ ~ 1.5 steradian). This latter constraint is consistent 
with the morphology of the outflow (e.g., |Goetz et al. 1990| ; [McKeith et al 19^ ; [Shop bell 



fc Bland-Hawthorn 1998]). Alternatively, we can use the empirical measurements of the 



gas pressure in the inner region of the M 82 nebula ( [Heckman, Armus, fc Miley 199C ) and 



extrapolate these outward, assuming that the wind ram pressure falls as r~^. The pressure 
at a radius of 1 kpc was measured to be ~ 5 x 10^^'^ dyne cm~^, and so the predicted wind 
ram pressure at a radius of 11.6 kpc would be ~ 4 x 10~^^ dyne cm~^, about a factor three 
smaller than the estimated value (for fx = 1). In view of the uncertainties in the estimated 
physical parameters for the X-ray ridge, we regard this level of disagreement as acceptable. 

If the optical line emission is also shock-excited by the outflowing wind, then the lu- 
minosity observed at Ha should be consistent with the wind energetics. We can estimate 
the amount of emission from a shock as follows. The total energy dissipated in a shock 
is, Lshock = l/2pv4wf = 1.0 X ergs s~\ where A^pc is the shock surface area 

in kpc^, uq is the preshock density (cm~^), and Vg is in units of km s~^. From [Binette 



Dopita, fc Tuohy (1985)| , we have using the equation above, L{Ha) shock ~ ^L{H(3)r ■ ~ 



shock 



3L,;,ocfc/100 = 1.7 yiW^nQAkpc{vs/S2Qf ergs s-\ if 90 km s'^ < < 1000 km s'^ Adopting 
the geometry above, the total surface area that the ridge would present to the wind would 
be 3.7 by 3.7 kpc. The observed Ha luminosity of the ridge (2.4 x 10^^ erg s~^) then re- 
quires that no(ws/820)'^ = 1.0 x 10"'^, while the pressure estimated in the X-ray ridge (1.3x 
10~^^ dynes cm~^) requires 7x0(^5/820)^ = 8 x 10""^. Thus, provided that the shock speed is 
greater than 100 km s~^, it can both produce the observed Ha emission-line luminosity and 
be consistent with the pressure estimated for the hot shocked gas observed in the X-rays. 



4.4. The Relationship Between the Ha and X-Ray Emission 



What is the physical relation between the Ha and X-ray emission in this type of model 
in which the superwind drives a shock into an ambient cloud in the halo of M 82? One 
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possibility is that the shock driven into the cloud by the wind is fast enough (i.e., ~ 800 km 
s~^) to produce the X-ray emission, and that the shocked gas then cools radiatively down 
to temperatures of-order 10^ K and produces the observed Ha emission via recombination. 
The difficulty with this simple model is that the estimated radiative cooling time in the 
X-ray gas (~ 3 x 10^ years) is very long compared to the characteristic dynamical time for 
the X-ray ridge: tdyn ~ (0.9 kpc/820 km s~^) ~ 1 x 10® years. Simply put, this model 
would require the X-ray ridge to be ~ 300 kpc thick and the Ha emission would be at the 
back edge of this very thick cooling zone! An additional persuasive argument against this 
model is provided by the rather quiescent kinematics of the Ha-emitting material reported 
by DB. Their optical spectra show that the Ha emission in the ridge spans a velocity range 
of only about 10^ km s~^, which is inconsistent with material cooling behind a shock with 
a velocity nearly an order-of- magnitude larger. While the relatively small difference DB 
report between the mean velocity of the ridge and M 82 proper (blueshift of 50 to 200 km 
s~^) could perhaps be explained as a projection effect in an outflow seen nearly in the plane 
of the sky, this effect could not explain the narrowness of the Ha emission-line in the ridge 
(~ 10^ km s~^) without appealing to an unreasonably idealized shock geometry (which is 
bellied by the morphological complexity of the Ha emission in the image published by DB). 

A much more plausible alternative is that we are actually observing two physically- 
related shocks in the M 82 ridge. We propose that as the superwind encounters a large cloud 
in the halo of M 82, it drives a relatively slow (of-order 10^ km s^^) radiative shock into 
the cloud. At the same time, this encounter also results in a fast (820 km s~^) stand-off 
bow-shock in the superwind fluid upstream from the cloud. The slow cloud shock (aided and 
abetted by ionizing radiation from M 82) produces the Ha emission. This is consistent with 
the quiescent Ha kinematics. The fast bow-shock in the wind produces the X-ray emission 
and is predicted to only be offset from front edge of the cloud by a fraction of the size of 
the cloud ( [Klein, McKee, fc Colella 1994| ; [Bedogni fc Woodward 1990| ). Pressure balance 



between the two shocks requires that PwindViow = PcioudVihock, clouds where pwind and pdoud are 
the preshock densities in the wind and cloud respectively. The observed X-ray pressure and 
temperature imply that ^bow ~ 820 km s~^ and n^ind^ion ~ 1 x 10^^ cm^^. Thus - for example 
- if the shock speed in the cloud is Vg^doud = 100 km s~^, then ndoud ~ 0.07 cm~^. We note 
further that the velocity of the bow shock is just the difference between the outflow velocity 
of the superwind itself and that of the shocked cloud. Assuming that latter velocity is much 
smaller than the former (as is suggested by the kinematic information provided by DB), the 
wind velocity implied by the X-ray temperature is ~ 820 km s~^. To drive a wind with a 
terminal velocity of 820 km s~^, hot gas would need to be injected into the starburst at a 
mean mass-weighted temperature of 1 x 10'^ K ( Chevalier fc Clegg 1985 ). This temperature 



is consistent with X-ray spectroscopy of the hot gas in the core of M 82 (DWH, [Moran & 



- 15 - 



Lehnert 1997; Strickland, Ponman, fc Stevens 1997 ). 



We have seen above that the radiative coohng time in the X-ray ridge is ~300 times 
longer than the characteristic dynamical time for the ridge. This implies that the wind bow 
shock is adiabatic. This can also be seen by comparing the total rate of energy dissipation in 
the shock {Lghock = ^ I '^PwindAridgevlow = 7x 10^° erg for the parameters estimated above) 
to the observed X-ray luminosity of the ridge (2.2 x 10'^^ erg s~^). Thus, the thickness of 
the X-ray ridge in the direction of the outflow is set by the dimensions of the region where 
adiabatic expansion and cooling causes the X-ray surface-brightness to drop precipitously. 
This region will be of-order the size of the obstacle in the flow that has caused the bow- 
shock to develop (the cloud). This is consistent with the observed X-ray morphology and the 
good spatial correlation seen between the Ha and X-ray emission. We have also estimated 
previously that the superwind can supply the ridge with its estimated thermal energy content 
in a timescale of-order 10® years. This is roughly equal to the dynamical time in the ridge 
(the timescale for adiabatic cooling), as required. 

We conclude that it is energetically feasible for the cloud to be shock-heated by a 
superwind driven by the starburst in M 82. A wind/cloud shock interpretation leads to no 
uncomfortable estimates of the luminosity, covering fraction, or time scales. In addition, as 
noted previously, there is ample evidence that M 82 lies in an "HI halo" ( |Yun et al. 1993 



1994). It is likely that the outflowing gas has intercepted and shock- heated some of this 
halo material, which is observed to be clumpy over large scales (see [Yun et al. 1993| ). We 



also note that perhaps one may have expected to find stronger evidence a priori for such 
interactions on the north side of the plane of M 82 compared to the south side: |Yun et aT 



1994| find that the velocity gradient in the HI is much higher on the northern side (12 km 
s~^ kpc^^) of the plane of M 82 compared to that on the southern side (3-6 km s~^ kpc^^). 
This higher gradient may be related to our observation of co-spatial Ha and thermal X-ray 
emission at large distances from the plane of M 82. 



4.5. Where is the Cloud? 

We have shown that a simple model of the interaction between a fast (~ 800 km s^^) 
wind and a massive cloud can quantitatively explain the physical and dynamical properties 
of the ridge. However, it a fair question to ask why the putative cloud is not observed in the 
sensitive HI map published by |Yun et al. 1993| . More specifically, we have estimated above 



that the high pressures in the X-ray gas and the low velocity of the shock driven into the 
cloud by the wind together mean that the preshock density of the cloud must be of-order 
10"^ cm~^. For a cloud with a line-of-sight dimension of 3.7 kpc, the implied cloud column 
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density is then ~10^^ cm ^. The HI map in |Yun et al. 1993 constrains the HI column density 
of the ridge to be < 2.7x 10^^ cm~^. This implies that the cloud must now be fully ionized. 

This could be accomplished by either the wind-driven shock or by the ionizing radiation 



cloud 



from M 82. The time for the shock to traverse the cloud will be ~ 10 years for Vshock 
100 km s~^, comparable to the estimated age of the starburst and its superwind. Similarly, 
the velocity of a photo-ionization front moving though the cloud is just Vi ~ ^Lyc^dlud (^-g-? 
|Oster brock 1989| ), where ^Lyc is the incident flux of ionizing radiation. The observed Ha 
luminosity and the assumed cross-sectional area the cloud presents to the starburst (3.7 by 
3.7 kpc) imply that the average value of ^Lyc is 1.5 xlO^/^J s~^ cm^^, where fats is the 
fraction of the incident ionizing photons absorbed by the cloud (i.e., allowing the cloud to 
be optically-thin). For a preshock density of 0.07 cm~'^ in the cloud, the implied velocity is 
Vi ~ 210 fabs~^ km s~^, and the time for the front to cross the cloud is about 4 x 10^ iabs~^ 
years. Thus it is plausible that an initially neutral cloud like those seen in the HI maps of 
Yun et al. (1993,1994) has indeed been fully ionized in the time since the starburst and its 
superwind turned on. 

Our rough estimates suggest that the timescale for the cloud to be fully photoionized is 
likely to be somewhat shorter than the timescale for it to be fully shock-ionized. If correct, 
this would mean that the Ha emission produced by the gas cooling and recombining behind 
the shock would have a high density and pressure, while the rest of the ionized cloud (the 
as-yet unshocked material) would be at a much lower pressure and density. We speculate 
that the bright knots and loops of Ha visible in the image in DB may correspond to shock- 
excited emission from gas with a high pressure, while the fainter and more diffuse emission 
may come from the low-pressure starburst-photoionized gas that has not yet been shocked. 



4.6. Implications for the "Superwind" Phenomenon 

It is often difficult to associate discrete X-rays sources seen in the halos of starburst 
galaxies with activity in the starburst galaxy itself. Often, it is assumed that most discrete 
objects observed in the halos of starbursts are more distant objects (usually QSOs) seen in 
projection. Finding a spatial association between Ha emission at the redshift of the starburst 
galaxy and X-ray emission is a irrefutable way of being able to associate the X-ray emission 
directly with activity within the starburst galaxy. 

Having argued that this co-spatial region of soft X-ray and Ha emission must be due to 
the outflowing superwind, it is now evident that at least some starburst galaxies are able to 
drive material to much larger distances from their nuclei than would be inferred using the 
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relatively high surface brightness X-ray or line emission that is generally observed along their 
minor axes (see e.g., |Lehnert fc Heckman 1996| ; |Heckman, Armus, fc Miley 1990| ; DWH). In 



M 82 for example, the inner "continuous" region of bright Ha and X-ray emission extends 
only about 5 arcminutes from the nucleus along the minor axis. The detached "ridge" 
of emission that we have discussed is at a projected angular separation from the nucleus 
of over twice this distance (about 11 arcminutes). Thus without making the association 
between the Ha and X-ray emission, a variety of interpretations for the source of the X- 
ray emission from this region would be possible. This obviously has important implications 
for our understanding of superwinds. It emphasizes that one must be careful in using the 
observed spatial scales of the extended emission in starbursts to either access whether or 
not the wind material is able to escape the potential of the host galaxy or to estimate the 
dynamical age of the outflow. 

It is interesting in this vein to consider the long-term fate of the hot gas in the X- 
ray ridge of M 82. To do so, we will compare the observed temperature of the gas to the 
"escape temperature" at that location in the M 82 halo. Following [Wang (1995)| , the escape 
temperature for hot gas in a galaxy potential with an escape velocity Vg^ is given by: T^s — 
1.1 X 10^ (ves/100 km s~^)^ K. The combined CO plus HI rotation curve of M 82 shows a 
roughly Keplerian decrease with radius for radii greater than 0.4 kpc, and Vrot has dropped 
to 60 km s~^ at a radius of 3.5 kpc ( |Sofue 1998 ). If this Keplerian fall-off continues to 



larger radii, the implied escape velocity at a radius of 11.6 kpc is only 45 km s~^ and the 
corresponding escape temperature is Tes ^ 2 x 10^ K. This is 450 times cooler than the 
observed X-ray emission. Even if we assume instead that M 82 has an isothermal dark 
matter halo with a depth corresponding to a circular velocity of 60 km s~^ (the maximum 
permitted by the directly measured rotation curve) and further assume that this halo extends 
to a radius of 100 kpc, the escape velocity at r = 11.6 kpc is about 150 km s~^ and Tes = 2.5 
xlO^ K. This is still a factor of 36 below the observed temperature of the gas in the ridge. 
This implies that the observed hot gas (and by inference, the superwind itself) will be able 
to escape the gravitational potential of M 82. 

The X-ray/Ha ridge has particularly interesting implications for the physics of super- 
wind emission. While the existence of superwinds is now well-established, the process(es) by 
which they produce X-ray emission is (are) not clear. This topic has been nicely reviewed 
by [Strickland 1998| . If the outflowing wind consists purely of the thermalized ejecta from 
supernovae and stellar winds, the wind fluid is hot (~ 10*^ K prior to adiabatic cooling) and 
is so tenuous that it will be an insignificant source of X-ray emission ( [Chevalier fc (Jlegg 



[1985[ ; [Suchkov et al. 1994[ ). Various alternatives have been proposed to boost the X-ray 
luminosity produced by a superwind: 1) Substantial quantities of material could be mixed 
into the wind in or near the starburst as the stellar ejecta interact with the starburst ISM 
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(the wind could be centrally "mass-loaded" - Suchkov et al. 1996 ; Hartquist et al. 1997 ). 2 



Large quantities of gas could be added in situ as the wind overtakes and then evaporates or 
shreds gas clouds in the galactic halo (e.g., [Suchkov et al. 19^1). These clouds could either 



be pre-existing clouds in the halo or material from the disk of the starburst galaxy that 
has been carried into the halo by the wind. 3) The tenuous wind fluid could drive a shock 
into a denser volume-filling galactic halo, with the observed X-ray emission arising from the 
shocked-halo material rather than the wind-fluid itself. 

In the case of the ridge of emission in the halo of M 82, we have argued that the 
observable X-ray emission arises as the wind fluid - which due to adiabatic expansion and 
cooling would otherwise have an X-ray surface-brightness below the level detectable in the 
existing ROSAT image - encounters a few-kpc-size-scale cloud in the halo. We have argued 
that the observed X-ray emission is produced in the bow-shock created in the superwind just 
upstream of the cloud. This process might be particularly relevant to the case of M 82, which 
is immersed in an extensive system of tidally-liberated gas clouds, shared with its interaction 
partners M81 and NGC3077. Such material would also provide a potential source of optical 
line emission by either intercepting a small portion of the ionizing radiation produced by 
the central starburst or as it is shock-heated by the superwind. Given the strong connection 
between galaxy interactions/mergers and the starburst phenomenon (e.g., [Sanders fc Mirabel 



1996 and references therein), this mechanism may be generally applicable to superwinds. 



5. Summary and Conclusions 

We report an analysis of a newly-discovered (pevine fc Bally 1999| ) region of spatially- 
coincident X-ray and Ha emission to the north of the prototypical starburst/superwind 
galaxy M 82. At the distance of M 82 (3.63 Mpc) the dimensions of this correlated ridge of 
emission are 3.7 x 0.9 kpc and it lies at a projected distance of about 11.6 kpc above the 
plane directly along the minor axis of M 82 (and hence along the axis of the superwind). 
We measure a total Ha flux from the ridge of 1.5 x 10^^^ ergs s^^ cm This flux yields 
a total Ha luminosity of 2.4 x 10'^^ ergs s~^. The total Ha flux for the whole of M 82 is 
approximately 4.6 x 10^^^ ergs s~^ cm Thus the Ha ridge constitutes approximately 
0.3% of the total observed Ha flux from the galaxy. X-ray emission is seen over the same 
region as the Ha emission by the ROSAT PSPC and HRI. With the 15 arcsec (260 pc) off- 
axis spatial resolution of the HRI, the ridge is diffuse and featureless except for a (probably 
unrelated) soft point source at its ENE end. The best flt to the X-ray emission from the 
ridge is a single thermal component with a temperature of kT=0.80±0.17 keV, absorbed by 
a column density = 3.7 x 10^° cm~^. The absorption is consistent with the foreground 
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Galactic column. The total unabsorbed flux from the ridge of emission is 1.4 x 10 ergs 
cm^^ s^^ and its himinosity of 2.2 x 10"^^ erg comprises about 0.7% of the total X-ray 
emission from M 82 in the 0.1 - 2.4 keV ROSAT band. 

We find that the observed Ha emission from the ridge could be produced via photoion- 
ization by the dilute radiation from the starburst propagating into the M 82 halo along the 
path cleared by the superwind. On the other hand, the X-ray emission from the ridge can 
be most naturally explained as the result of shock-heating associated with an encounter be- 
tween the starburst-driven galactic "superwind" and a large photoionized cloud in the halo 
of M 82. Shock-heating may also contribute significantly to the excitation of the Ha emis- 
sion, especially the regions of the highest surface brightness. We find that the Ha and X-ray 
luminosities and the X-ray temperature and estimated gas pressure can all be quantitatively 
understood in this context. The specific model we favor is one in which the encounter be- 
tween the superwind and the cloud drives a slow radiative shock into the cloud (producing 
some of the Ha emission) and also leads to an adiabatic X-ray-emitting bow-shock in the 
superwind fluid upstream from the cloud. The cloud is now fully ionized, and so is not 
detected in the HI maps of Yun et al. (1993, 1994). The wind-cloud interaction process may 
be especially relevant to M 82, whose halo contains tidally-liberated gas clouds. However, 
it may well occur generally in powerful starbursts (which are usually triggered by galaxy 
interactions or mergers). 

The existence of a region of spatially coincident X-ray and Ha emission with strong 
evidence for wind/cloud interaction is important for several reasons. Firstly, our analysis 
sheds new light on the physical processes by which galactic winds produce X-ray emission, 
and implies that wind/cloud collisions are an important emission-producing process. The 
emission ridge in M 82 highhghts the potential pitfalls in using the observed X-ray or Ha 
emission to trace the size and infer the ages or fates of galactic winds: the winds may extend 
to radii far larger than the region over which they have densities and temperatures high 
enough to produce observable emission. Wind-cloud encounters then offer us the possibility 
of 'lighting up' the wind at radii that are otherwise inaccessible. The X-ray ridge in M 82 
implies that the wind is able to reach large heights above the plane of the galaxy, thereby 
increasing the likelihood that it will be able to eventually escape the galaxy potential well. 
Indeed, the observed gas temperature exceeds the local "escape temperature" from the M 82 
gravitational potential well by about two orders-of-magnitude. This result strengthens the 
case that starburst-driven galactic winds are important sources of the metal-enrichment and 
heating of the Intergalactic Medium. 
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FIGURE CAPTIONS 

Fig. 1. — The greyscale is the Ha image and the contours are the ROSAT PSPC image of 
M 82. The lowest surface brightness emission visible in the Ha corresponds to a flux of 3.5 x 
10~^^ ergs s~^ cm ~^ arcsec"^. The spatial coincidence of the Ha and X-ray emission is quite 
good. The images were aligned using using the alignment determined for the continuum 
image. The coordinates are shown for J2000. 

Fig. 2. — The greyscale is a narrow-band continuum image (central wavelength approximate 
6658A) image and the contours are the smoothed {a=10" gaussian) ROSAT PSPC image of 
M 82. The images were ahgned using what is hkely to be a background AGN and galaxy that 
were both visible in the optical continuum and X-ray images. The coordinates are shown for 
J2000. 

Fig. 3. — The Rosat PSPC spectrum of the diffuse X-ray emission from the ridge to the right 
of the point source (circles) and superimposed on the point source (crosses). Data below 0.6 
keV are ignored for the latter case to eliminate the point source. Also shown are the best-fit 
MEKAL plasma models (see text) folded through the instrumental response. 
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